Group Analysis of Self-organizing Maps based on Functional MRI using Restricted Frechet Means

Studies of functional MRI data are increasingly concerned with the estimation of differences in spatio-temporal networks across groups of subjects or experimental conditions. Unsupervised clustering and independent component analysis (ICA) have been used to identify such spatio-temporal networks. While these approaches have been useful for estimating these networks at the subject-level, comparisons over groups or experimental conditions require further methodological development. In this paper, we tackle this problem by showing how self-organizing maps (SOMs) can be compared within a Frechean inferential framework. Here, we summarize the mean SOM in each group as a Frechet mean with respect to a metric on the space of SOMs. We consider the use of different metrics, and introduce two extensions of the classical sum of minimum distance (SMD) between two SOMs, which take into account the spatio-temporal pattern of the fMRI data. The validity of these methods is illustrated on synthetic data. Through these simulations, we show that the three metrics of interest behave as expected, in the sense that the ones capturing temporal, spatial and spatio-temporal aspects of the SOMs are more likely to reach significance under simulated scenarios characterized by temporal, spatial and spatio-temporal differences, respectively. In addition, a re-analysis of a classical experiment on visually-triggered emotions demonstrates the usefulness of this methodology. In this study, the multivariate functional patterns typical of the subjects exposed to pleasant and unpleasant stimuli are found to be more similar than the ones of the subjects exposed to emotionally neutral stimuli. Taken together, these results indicate that our proposed methods can cast new light on existing data by adopting a global analytical perspective on functional MRI paradigms.Comment: 23 pages, 5 figures, 4 tables. Submitted to Neuroimag

The Spatial Extent and Distribution of Star Formation in 3D-HST Mergers at z~1.5

We present an analysis of the spatial distribution of star formation in a sample of 60 visually identified galaxy merger candidates at z>1. Our sample, drawn from the 3D-HST survey, is flux-limited and was selected to have high star formation rates based on fits of their broad-band, low spatial resolution spectral energy distributions. It includes plausible pre-merger (close pairs) and post-merger (single objects with tidal features) systems, with total stellar masses and star formation rates derived from multi-wavelength photometry. Here we use near-infrared slitless spectra from 3D-HST which produce Halpha or [OIII] emission line maps as proxies for star-formation maps. This provides a first comprehensive high-resolution, empirical picture of where star formation occurred in galaxy mergers at the epoch of peak cosmic star formation rate. We find that detectable star formation can occur in one or both galaxy centres, or in tidal tails. The most common case (58%) is that star formation is largely concentrated in a single, compact region, coincident with the centre of (one of) the merger components. No correlations between star formation morphology and redshift, total stellar mass, or star formation rate are found. A restricted set of hydrodynamical merger simulations between similarly massive and gas-rich objects implies that star formation should be detectable in both merger components, when the gas fractions of the individual components are the same. This suggests that z~1.5 mergers typically occur between galaxies whose gas fractions, masses, and/or star formation rates are distinctly different from one another.Comment: Accepted for publication in MNRAS, 16 pages, 10 figure

No More Active Galactic Nuclei in Clumpy Disks Than in Smooth Galaxies at z~2 in CANDELS / 3D-HST

We use CANDELS imaging, 3D-HST spectroscopy, and Chandra X-ray data to investigate if active galactic nuclei (AGNs) are preferentially fueled by violent disk instabilities funneling gas into galaxy centers at 1.3<z<2.4. We select galaxies undergoing gravitational instabilities using the number of clumps and degree of patchiness as proxies. The CANDELS visual classification system is used to identify 44 clumpy disk galaxies, along with mass-matched comparison samples of smooth and intermediate morphology galaxies. We note that, despite being being mass-matched and having similar star formation rates, the smoother galaxies tend to be smaller disks with more prominent bulges compared to the clumpy galaxies. The lack of smooth extended disks is probably a general feature of the z~2 galaxy population, and means we cannot directly compare with the clumpy and smooth extended disks observed at lower redshift. We find that z~2 clumpy galaxies have slightly enhanced AGN fractions selected by integrated line ratios (in the mass-excitation method), but the spatially resolved line ratios indicate this is likely due to extended phenomena rather than nuclear AGNs. Meanwhile the X-ray data show that clumpy, smooth, and intermediate galaxies have nearly indistinguishable AGN fractions derived from both individual detections and stacked non-detections. The data demonstrate that AGN fueling modes at z~1.85 - whether violent disk instabilities or secular processes - are as efficient in smooth galaxies as they are in clumpy galaxies.Comment: ApJ accepted. 17 pages, 17 figure

Differential Neural Responses to Food Images in Women with Bulimia versus Anorexia Nervosa

BACKGROUND: Previous fMRI studies show that women with eating disorders (ED) have differential neural activation to viewing food images. However, despite clinical differences in their responses to food, differential neural activation to thinking about eating food, between women with anorexia nervosa (AN) and bulimia nervosa (BN) is not known. METHODS: We compare 50 women (8 with BN, 18 with AN and 24 age-matched healthy controls [HC]) while they view food images during functional Magnetic Resonance Imaging (fMRI). RESULTS: In response to food (vs non-food) images, women with BN showed greater neural activation in the visual cortex, right dorsolateral prefrontal cortex, right insular cortex and precentral gyrus, women with AN showed greater activation in the right dorsolateral prefrontal cortex, cerebellum and right precuneus. HC women activated the cerebellum, right insular cortex, right medial temporal lobe and left caudate. Direct comparisons revealed that compared to HC, the BN group showed relative deactivation in the bilateral superior temporal gyrus/insula, and visual cortex, and compared to AN had relative deactivation in the parietal lobe and dorsal posterior cingulate cortex, but greater activation in the caudate, superior temporal gyrus, right insula and supplementary motor area. CONCLUSIONS: Women with AN and BN activate top-down cognitive control in response to food images, yet women with BN have increased activation in reward and somatosensory regions, which might impinge on cognitive control over food consumption and binge eating

REQUIEM-2D: Spatially Resolved Stellar Populations from HST 2D Grism Spectroscopy

We present a novel Bayesian methodology to jointly model photometry and deep Hubble Space Telescope (HST) 2d grism spectroscopy of high-redshift galaxies. Our requiem2d code measures both unresolved and resolved stellar populations, ages, and star-formation histories (SFHs) for the ongoing REQIUEM (REsolving QUIEscent Magnified) Galaxies Survey, which targets strong gravitationally lensed quiescent galaxies at z~2. We test the accuracy of \texttt{requiem2d} using a simulated sample of massive galaxies at z~2 from the Illustris cosmological simulation and find we recover the general trends in SFH and median stellar ages. We further present a pilot study for the REQUIEM Galaxies Survey: MRG-S0851, a quintuply-imaged, massive ($\log M_* / M_\odot = 11.02 \pm 0.04$) red galaxy at $z=1.883\pm 0.001$. With an estimated gravitational magnification of $\mu = 5.7^{+0.4}_{-0.2}$, we sample the stellar populations on 0.6 kpc physical size bins. The global mass-weighted median age is constrained to be $1.8_{-0.2}^{+0.3}$ Gyr, and our spatially resolved analysis reveals that MRG-S0851 has a flat age gradient in the inner 3 kpc core after taking into account the subtle effects of dust and metallicity on age measurements, favoring an early formation scenario. The analysis for the full REQUIEM-2D sample will be presented in a forthcoming paper with a beta-release of the requiem2d code.Comment: 29 pages, 19 figures, accepted for publication in Ap

The DEEP2 Galaxy Redshift Survey: Design, Observations, Data Reduction, and Redshifts

We describe the design and data sample from the DEEP2 Galaxy Redshift Survey, the densest and largest precision-redshift survey of galaxies at z ~ 1 completed to date. The survey has conducted a comprehensive census of massive galaxies, their properties, environments, and large-scale structure down to absolute magnitude M_B = -20 at z ~ 1 via ~90 nights of observation on the DEIMOS spectrograph at Keck Observatory. DEEP2 covers an area of 2.8 deg^2 divided into four separate fields, observed to a limiting apparent magnitude of R_AB=24.1. Objects with z < 0.7 are rejected based on BRI photometry in three of the four DEEP2 fields, allowing galaxies with z > 0.7 to be targeted ~2.5 times more efficiently than in a purely magnitude-limited sample. Approximately sixty percent of eligible targets are chosen for spectroscopy, yielding nearly 53,000 spectra and more than 38,000 reliable redshift measurements. Most of the targets which fail to yield secure redshifts are blue objects that lie beyond z ~ 1.45. The DEIMOS 1200-line/mm grating used for the survey delivers high spectral resolution (R~6000), accurate and secure redshifts, and unique internal kinematic information. Extensive ancillary data are available in the DEEP2 fields, particularly in the Extended Groth Strip, which has evolved into one of the richest multiwavelength regions on the sky. DEEP2 surpasses other deep precision-redshift surveys at z ~ 1 in terms of galaxy numbers, redshift accuracy, sample number density, and amount of spectral information. We also provide an overview of the scientific highlights of the DEEP2 survey thus far. This paper is intended as a handbook for users of the DEEP2 Data Release 4, which includes all DEEP2 spectra and redshifts, as well as for the publicly-available DEEP2 DEIMOS data reduction pipelines. [Abridged]Comment: submitted to ApJS; data products available for download at http://deep.berkeley.edu/DR4

The 3D-HST Survey: <i>Hubble Space Telescope</i> WFC3/G141 Grism Spectra, Redshifts, and Emission Line Measurements for ~ 100,000 Galaxies

We present reduced data and data products from the 3D-HST survey, a 248-orbit $HST$ Treasury program. The survey obtained WFC3 G141 grism spectroscopy in four of the five CANDELS fields: AEGIS, COSMOS, GOODS-S, and UDS, along with WFC3 $H_{140}$ imaging, parallel ACS G800L spectroscopy, and parallel $I_{814}$ imaging. In a previous paper, we presented photometric catalogs in these four fields and in GOODS-N, the fifth CANDELS field. Here we describe and present the WFC3 G141 spectroscopic data, again augmented with data from GO-1600 in GOODS-N (PI: B. Weiner). We developed software to automatically and optimally extract interlaced two-dimensional (2D) and one-dimensional (1D) spectra for all objects in the Skelton et al. (2014) photometric catalogs. The 2D spectra and the multi-band photometry were fit simultaneously to determine redshifts and emission line strengths, taking the morphology of the galaxies explicitly into account. The resulting catalog has redshifts and line strengths (where available) for 22,548 unique objects down to ${{JH}}_{\mathrm{IR}}\leq 24$ (79,609 unique objects down to ${{JH}}_{\mathrm{IR}}\leq 26$). Of these, 5459 galaxies are at $z > 1.5$ and 9621 are at $0.7< z< 1.5$, where Hα falls in the G141 wavelength coverage. The typical redshift error for ${{JH}}_{\mathrm{IR}}\leq 24$ galaxies is ${\sigma }_{z}\approx 0.003\times (1+z)$, i.e., one native WFC3 pixel. The $3\sigma$ limit for emission line fluxes of point sources is $2.1\times {10}^{-17}$ erg $s^{-1} cm^{-2}$. All 2D and 1D spectra, as well as redshifts, line fluxes, and other derived parameters, are publicly available